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Abstract 

We evaluate differential distributions for the four-body pp — t- nmr'^TT^ reaction. The ampli- 
tude for the process is calculated in the Regge approach including many diagrams. We make 
predictions for possible future experiments at RHIC and LHC energies. Very large cross sec- 
f^ ' tions are found which is partially due to interference of a few mechanisms. Presence of several 

^^ ' interfering mechanisms precludes extraction of the elastic ir'^ir^ scattering cross section. Absorp- 

S ■ tion effects are estimated. Differential distributions in pseudorapidity, rapidity, invariant two-pion 

mass, transverse-momentum and energy distributions of neutrons are presented for proton-proton 
, collisions at ^ = 500 GeV (RHIC) and ^/s = 0.9, 2.36 and 7 TeV (LHC). Cross sections with 

^ I experimental cuts are presented. 

^; PACS numbers: 11.55.Jy, 13.85.Lg, 14.20.Dh 
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I. INTRODUCTION 

The total and elastic cross sections are basic objects of the scattering theory. While 
the proton-proton, proton-antiproton or pion-proton can be directly measured (see e.g.|l|) 
the pion-pion scattering is not directly accessible. It was suggested recently [3] how to 
extract the total tt+tt^ cross section in the high-energy region. Here it was suggested to use 
scattering of virtual tt+'s which couple to the nucleons with well known coupling constant 
and are subsequently promoted by the interaction onto their mass shell in the final state. 
The final pions are then associated with outgoing neutrons. 

Can a similar method be used to extract the elastic tt~^tt~^ scattering by analysis of the 
pp — )■ nrnr^'K^ reaction? We wish to address this issue in the present paper. ^ The energy 
dependence of the total and possibly elastic cross section of pion-pion scattering would 
be very useful and supplementary information for the groups which model hadron-hadron 
interactions in the soft sector (see e.g. [^). 

It was realized over the last decade that the measurement of forward particles can be 
an interesting and useful supplement to the central multipurpose LHC detectors (ATLAS, 
CMS). The main effort concentrated on the design and construction of forward proton 
detectors [q. Also Zero-Degree Calorimeters (ZDC's) haA/'e been considered as a useful 
supplement. It will measure very forward neutrons and photons in the pseudorapidity region 
\7]\ > 8.5 at the CMS (gJ (see also [7(]) and the ATLAS ZDC's provide coverage of the 
region \r]\ > 8.3 |8|. It was shown recently that the CMS (Compact Muon Spectrometer) 
Collaboration ZDC's provide a unique possibility to measure the 7r+7r+ total cross section 

0. 

Even at high-energy the major part of the phase space of a few-body reactions is populated 
in soft processes which cannot be calculated within perturbative QCD. Only limited corner of 
the phase space, where particles are produced at large transverse momenta, can be addressed 
in the framework of pQCD. At high energy the Regge approach is the most efficient tool 
to describe total cross section, elastic scattering as well as different 2 — ;■ 2 reactions |9|. In 
the present paper we shall show how to construct the amplitude for the considered 2 — )■ 4 
process in terms of several 2 — )■ 2 soft amplitudes known from the literature. In the present 
analysis we will also include absorption effects as was done recently for three-body processes 
0. 

In the present paper we consider an example of an exclusive reaction with two forward 
neutrons. Given the experimental infrastructure the pp — )■ nn7r~^7r~^ is one of the reactions 
with four particles in the final state which could be addressed at LHC. 

II. AMPLITUDE OF EXCLUSIVE pp -^ nmT+TT+ REACTION 

A. Dominant difFractive amplitude 

The diffractive mechanisms involving pomeron and reggeon exchanges included in the 
present paper are shown in Fig{T] (with the four- momenta Pa + Pb -^ Pi + P2 + P3 + Pi)- 



^ After first version of our paper had been completed, a paper has appeared which also discusses the 
possibility of extraction of elastic 7r+7r+ cross section 3]. In our analysis we take into account many more 
possible mechanisms for the pp 



In principle in all diagrams shown the intermediate nucleon can be replaced by nucleon 
excited states. It is known that diffractive excitation of nucleons to inelastic states is rather 
large and constitutes about 1/3 of the elastic scattering. This number is, however, not 



relevant in our case, as it is to large extend due to the Deck type mechanism [ll[ which is 
included explicitly in our calculation. The remaining excitation to discrete nucleon states is 
rather small and difficult to calculate. A microscopic calculation must unavoidably include 
not only the structure of the nucleon but also of the nucleon excited states. The cross 
section for pp — )■ p + Nti-k of our interest is, however, only a fraction of mb [12|. That the 
contribution of excited discrete state is small can be also seen in the following way. First 
of all the diffractive transitions to discrete excited states are known to be much weaker 
than the elastic one. Secondly the gNN*-K coupling constants are much smaller than the 
dNN-K coupling constant [l3|. Finally the exact strength of the diffractive transitions are 
not known phenomenologically. Therefore in the following we neglect the contributions of 
diagrams with excited nucleon states. 
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FIG. 1: Diagrams for the exclusive production of vr^vr^ in -pp collisions at high energies. The 
stars attached to vr"^, n and p denote the fact they are off-mass-shell. ka~hi are four- vectors of the 
exchanged pomerons. 

Similarly as for the pp — > A^iV/o(1500) [iJ] and "pip — )■ ppTr+vr^ [jtp — )■ ppTr+vr^) 15|, [l6| 
reactions the amplitudes can be written in terms of pomeron (reggeon)-exchanges. Then the 
amplitude squared, averaged over the initial and summed over the final polarization states, 
for the 'p'p — )■ nm^^'n^ process can be written as: 

^w = ^ E iM:UA.A.+-+Ml-.Aa.r (2-1) 

AaAi,AlA2 

It is straightforward to evaluate the contribution shown in FiglH The diagrams will be 
called a) - i) for brevity. If we assume the 175 type coupling of the pion to the nucleon then 



the Born amplitudes read: ^ 

■^\!\,^\l\2 = u{Pl^^l)h5Sn{tl)u{pa,Xa)V2gnNNF^NNiti) 

X u{p2,\2)h5Sn{t2)u{pb,\b)V2g^NNF7TNN{t2) , (2.2) 



(b) 



-^J^Afc^AiAa = M(Pl,Ai)i75S'^(ti)M(p„, Aa)y25f^ArArF^7VAr(ti) 



X F;//(tO i..3C- (^j (^] exp [^^klj F:ffiu2) 

X M(p2,A2)i755'n(M2)M(Pb, Afe)y25f^ArArF^ArAr(M2) , (2.3) 



-^Sa6^AiA2 = u{pi,Xi)i-f5Sn{Ui)u{pa,Xa)V2g^NNF^NN{Ul) 

/ \aip(fc?)-l / \ aiv(Mi)-| /n \ 

X F^ffiu,) i.MC- (^j (|ij exp [^^kVj F;//(t,) 

X u{p2,\2)h5S7r(t2)u{pb,\b)V2g7rNNF^NN(t2) , (2.4) 

•^AaAb^AiAa = M(pi,Ai)i75S'„(Mi)M(pa, Aa)v^fi'7rAfAr-F7r7VAf(Ml) 



X F^f^iui) iSi2Cjp 



I „ \ "ip(^d)-l /„ \ Q!iv(«i)-| /„ \ Q!jv(m2)-| 
NN / -^12 \ / -^IS \ / "^24 ^ 



So / \ "5t/i / V St/i 



X exp (^A:^) F://(.,) 

X M(j92,A2)i75'S'„(M2)M(p6, Ab)V25f^Ar7vF^Ar7v(M2) , (2.5) 

A^aJa,^AiA2 = w(Pl>Al)i755'^(ti)M(p„, Aa)y2fi',r7V7vF^ArAr(tl) 

X F;^^(ti)is234C^^(^) "''''" exp (^A:,^) F;//(s24) 

X M(p2,A2)i755'p(s24)M(Pfe, Afe)y25'^Ar7vF^ArAr(s24) , (2.6) 

^\a\b^\l\2 = M(Pl,Al)i755'„(Mi)M(pa,Aa)y2fi',rArAf-^7rAfAf(Ml) 

X F://(.0 i.,24Cr (^) (^) exp (^A:Jj F;//(s24) 

X M(p2,A2)i755'p(s24)M(P6, Afe)y25f^ArAfF^ArAr(s24) , (2.7) 



^ We show explicitly only amplitudes for pomeron exchange. The amplitudes for reggeon exchange can be 
obtained from those for pomeron exchange by replacing propagators by signature factors and trajectories. 



•^\!x,-,\^\2 = u{puXi)h5Sp{Si3)u{pa,Xa)V2g^NNF^NN{si3) 

X F;//(.,3) is,uC-/ (^ j exp [^^e^j F:ff{t,) 

X u{p2,X2)h5Sn(t2)u{pb,Xb)V2g^NNF^NN{t2), (2.8) 

-^£a6^AiA2 = M(pi,Ai)i75S'p(Si3)M(pa,Aa)y25'7r7V7V^7r7V7v(Sl3) 

X F;"(..3) i.,.3Cr (^) (^) exp (^*j) f.:"(,.,) 

X u{P2,X2)h5Sn{u2)u{pb,Xb)V2gTrNNF^NN{u2) , (2.9) 

-^£a6^AiA2 = '«(Pl5-^l)i75'S'p(Si3)u(p„,Aa)y25f^Ar7vF^ArAr(si3) 

X M(p2,A2)i75'S'p(s24)M(Pfe, Ab)y25f^7VArF^ArAr(s24) , (2.10) 

where the energy scale Sq is fixed at Sq = 1 GeV^ and Sth = ijn^ + m^)^. 

In the above equations u{pi,Xi), u{pf,Xf) = u\pf,Xf)'^^ are the Dirac spinors (nor- 
mahzed as u{p)u{p) = 2mAr) of the initial protons and outgoing neutrons with the four- 
momentum p and the helicities of the nucleons A. The propagators of virtual particles can 
be written as 

S„(.,,) = ""■■^-^^ + f"' . (2.12) 

s^is,) ^ i(!M.2:±^ , (2.13) 

where ti^2 = {Pa,b — Pi,2)^ and Ui^2 = {Pa,b — ^3,4)^ = "^12 ^^^ the four-momenta squared of 
transferred pions and neutrons, respectively ^. Sij = {pi+PjY = sf^ are the squared invariant 
masses of the (i, j) system, m^ and m„, rup are the pion and nucleons masses, respectively. 
The factor q-^nn is the pion nucleon coupling constant which is relatively well known [l7| 
((?^^jy/47r = 13.5 - 14.6). In our calculations the coupling constant is taken as g^j^^/An = 
13.5. 

Using the known strength parameters for the NN and ttN scattering fitted to the corre- 
sponding total cross sections (the Donnachie-Landshoff model 18|) we obtain C^^ , C^ and 



assuming Regge factorization [19|] Cp. The pomeron/reggeon trajectories determined from 
elastic and total cross sections are given in the linear approximation ^ (aj(t) = aj(0) + a^t) 
where the values of relevant parameters (the intercept aj(0) and the slope of trajectory a^ in 



^ In the following for brevity we shall use notation ti.2 which means ii or ^2- 

^ For simplicity we use the linear pomeron/reggeons trajectories, but further improvements are possible. 



TABLE I: Parameters of reggeon exchanges used in the present calculations. 
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GeV~^) are also taken from the Donnachie-Landshoff model 18| for consistency. Parameters 
of reggeon exchanges used in the present calculations are listed in Table [11 
The slope parameter can be writen as 



B{s) = Bo + 2ajp\n[ — 



(2.14) 



where Bq is the t-slope of the elastic differential cross section. In our calculation we use Bq: 
ByrN = 6.5 GeV~^, Bnn = 9 GeV~^ and Bj,^ = 4 GeV~^. The value of B.^^, is not well known, 



2i?^Ar — Bnn uM ■ We have parametrized the 



however the Regge factorization entails -B,, 

kl, ...,k'^ dependences in the exponential form (see formulas (12. 2p - (I2.10p ). 

We improve the parametrization of the amplitudes for neutron exchange (12. 3[ 

ajv("l,2)- 



EH IMD by the factors ^ 



Sth 



to reproduce the high-energy Regge dependence. The 



-2 



degenerate nucleon trajectory is aAr(Mi,2) = —0.3 + Q;^ni,2, with a'^- = 0.9 GeV 

The extra correction factors F!^^fj{k'^) (where fc^ = ti^2, mi,2, Sij) are due to off-shellness of 
particles. In the case of our 4-body reaction rather large transferred four-momenta squared 
k"^ are involved and one has to include non-point-like and off-shellness nature of the particles 
involved in corresponding vertices. This is incorporated via Fj^NNik'^) vertex form factors. 
We parametrize these form factors in the following exponential form: 



Fiti,2) 

^(Wl,2) 
F(S:,) 
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ti,2 - m; 

A2 



Mi,2 - m^ 

exp I 7^ I , 



A2 



exp 



[s,j - m;j 



A2 



(2.15) 
(2.16) 
(2.17) 



While four-momenta squared of transferred pions ti,2 < 0, it is not the case for transferred 
neutrons where Ui^2 < i^n- ■'■^ general, the cut-off parameter Aoff is not known but in 
principle could be fitted to the (normalized) experimental data. From our general experience 
in hadronic physics we expect Aoff ~ 1 GeV. Typical values of the ttNN form factor 



parameters used in the meson exchange models are A = 1.2-1 A GeV 20(|, however the 
Gottfried Sum Rule violation prefers smaller A ^ 0.8 GeV [21!]. In our calculation, if not 
otherwise mentioned, we use A = Aq// = 1 GeV. We shall discuss how uncertainties of the 
form factors influence our final results. 



B. Single and double charge exchanges with subleading reggeons p"*", a^ 

We wish to include also specific processes with isovector reggeon exchanges. We in- 
clude processes shown in FigjJl These processes involve p^ p^ — t- 7r+7r+ and a^a^ — )■ 7r+7r+ 
subprocesses. Unfortunately these subprocesses (or the reverse ones) could not be studied 
experimentally. 

The relevant coupling constants in diagrams b) and c) are not known and cannot be 
obtained from first principles and one has to refer to other reactions involving the same 
coupling constants. Such reactions are e.g. n^p — )■ Og P (where both P ^ p'^ exchanges are 
possible), Ti~p — ;■ a^^, Ti~p — )■ uj^n (only p^-reggeon exchange come into game), ir^p — )• p^p 
(7r°, uj^- and a^-reggeon exchanges) and 7r~p — )■ p^n (tt^, a^-reggeon exchanges). The details 
how to fix parameters of these two-body reactions are described in the Appendix. 




/ n(Pi) 





pW_^ 


7 




■+(Pa) 


' n 


"^-J 


.-ir+(pj 


■+(p..) 


-^ rT 


rl""^"" 


■-T+fe 




pil't) "^ 


^ 




c) 




V nlp2) 





FIG. 2: Diagrams with subleading charged reggeon exchanges in pp colHsions at high energies. 

The diagram a) in Fig|2] is topologically identical to the dominant diagram for the 
pp — )■ ppn^Ti^ reaction [16]. There, however, the pomeron-pomeron, pomeron-reggeon and 
reggeon-pomeron exchanges are the dominant processes. In addition to diagram a) there is 
possible also another mechanism with the intermediate pion replaced by a virtual photon. 
Because it requires two electromagnetic couplings instead of two strong couplings its contri- 
bution should be small. Because of the extra photon propagator it could be enhanced when 
H — )■ 0. However then the vertices should tend to zero. Therefore we can safely omit such 
a diagram. 

We write the amplitudes for the diagrams in Figl2] as: 

(, \ |Ai-Aa|/2 / \ aiR(ti)-l /R \ 

X ^(S34,ta) 

-|- crossedterm , (2.18) 

where ^(^34, to) refers to the central part of the diagrams 

1 



X 



A"^'is,,, Q = F:ff{Q ,/C- -— -^ ,/Cr F:ff{Q , (2.19) 



Areggeon-exch. ^ ^ j. \ V ^^ „ „ / „J \2 / '^34 \ ' ° (BmM_^ \ V 4 



^reggeon-e.c..^^^^^ t.) = ^ ^. ^34 (^}^.)^ ( ^ 1 exp ( "f^t J ^ .(2.20) 



TABLE II: Different realizations of diagram a) in FigJSj 
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In actual calculations we take Bmn = -Bttat and Bmm = -Btttt- Since, in the diagrams in 
Fig|2] and FigJH] we have reggeon exchanges rather than meson exchanges therefore formulas 
fl2.18t I2.20p give rather upper limit for the cross section. 

The parameterization of the amplitudes with subleading charged reggeon exchanges can- 
not be used in the region of resonances in ttN or/and vrvr subsystems [l6| . Therefore, the 
amplitude used in the calculations must contain restrictions on the four-body phase space. 
To exclude the regions of resonances we modify the parameterization of the amplitudes 
fl2.18p by multiplying cross section by a purely phenomenological smooth cut-off correction 
factor (see |l6|): 



Jcont [ ^^■kN/h 



exp (- 



< W-Wo \ 

1 + exp (M))" 



(2.21) 



The parameter Wq gives the position of the cut and the parameter a describes how sharp 
is the cut off. For large energies fcont''^{^-KN/TT-K) ~ 1 and close to kinematical threshold 

/cont '"(^TrAT/TTTr) ~ 0. In our calculatiou we take Wq = 2 GeV and a = 0.2 GeV. 

There is another class of diagrams shown in Figl3l The diagram (a) represents a generic 
amplitude, with particle sets (A, B, C) collected in Table HIl In contrast to the diagrams 
shown in Figj2] here both pions and subleading reggeons couple to nucleons. We shall not 
present explicit formulae for the corresponding amplitudes here. We shall show separate 
contributions of those processes in the Result Section. 
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FIG. 3: Diagrams with subleading reggeon at exchange in pp collisions at high energies. 
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FIG. 4: Schematic diagram for absorption effects due to proton-proton interaction. 
C. Absorptive corrections 



The absorptive correction in FigJH are calculated as described in [10| for the three body 
processes. Here the absorptive correction to the bare amplitude (see FigH]) can be writen 
as: 

SM,^x,^,,,Mt,P2t) = i/0^^A^ll^A,A.(P:* -Pi*>p7* -P2*) , (2.22) 

where Pa = Pa — h, Pb = Pb + h with momentum transfer kt- Above A^aJa -J.A1A2 ^^ ^ 
bare amplitude calculated as described in the previous subsections. T(s, kf) is an elastic 
proton-proton amplitude for the appropriate energy. It can be conveniently parametrized 

as: 

T{s, kl) = A^{s) eM-BNNkl/2) . (2.23) 

From the optical theorem we have Imy4o(s) = scrffj(s) (the real part is small in the high 
energy limit). Again the Donnachie-Landshoff parametrization [18| of the total pp or pp 
cross sections can be used to calculate the rescattering amplitude. 

In our analysis the vr+n interactions are not taken into account. They would further 
decrease the cross section. Given other theoretical uncertainties (form factors) it seems not 
worthy to take over the effort of performing very time-consuming calculations. Absorption 
effects for exclusive Higgs production are discussed e.g. in Ref . |24| . 

The cross section is obtained by assuming a general 2 — )■ 4 reaction: 

a = j -\M\^{2n) 5 (P^ + P^ " Pi - P2 - Ps -P4) ^^^^32j,, ^2.^32,,, (2^)32^,3 (2vr)32i?, - 

(2.24) 

To calculate the total cross section one has to calculate 8-dimensional integral numerically. 
The details how to conveniently reduce the number of kinematical integration variables are 
given elsewhere [16[. 

III. RESULTS 

We shall show our predictions for the pp — )■ nn7i~^iT~^ reaction for several differential 
distributions in different variables at selected center-of-mass energies W = 500 GeV (RHIC) 



TABLE III: Full-phase-space integrated cross section (in nib) for exclusive nmr^ir^ production at 
selected center-of-mass energies and different values of the form factor parameters. In parentheses 
we show cross sections including absorption effects. 



W = 0.5 TeV 



W = 0.9 TeV 



W = 2.36 TeV 



W = 7 TeV 



A = 0.8 GeV, Aoff = 1 GeV 

A = Aoff = 1 GeV 
A = 1.2 GeV, Aoff = 1 GeV 



0.34 (0.15) 
0.84 (0.37) 
1.45 (0.62) 



0.38 (0.16) 
0.95 (0.39) 
1.64 (0.66) 



0.47 (0.18) 
1.16 (0.42) 
2.01 (0.71) 



0.59 (0.19) 
1.47 (0.46) 
2.55 (0.77) 



and W = 0.9, 2.36 and 7 TeV (LHC). The cross section slowly rises with incident energy. In 
general, the higher energy the higher absorption effects. The results depend on the value of 
the nonperturbative, a priori unknown parameter of the form factor responsible for off-shell 
effects. In Table [TTll we have collected integrated cross sections for selected energies and 
different values of the model parameters. We show how the uncertainties of the form factor 
parameters affect our final results. 

In Figj5]we show distributions in pseudorapidity (r] = — ln(tan|), where 9 is the angle 
between the particle momentum and the beam axis) for the pp — >■ nn7r~^7r~^ reaction. The 
discussed reaction is very unique because not only neutrons but also pions are produced 
dominantly in very forward or very background directions forming a large size gap in pseu- 
dorapidity between the produced pions, about 12 units at W = 7 TeV. While neutrons can 
be measured by the ZDC's the measurement of very forward/backward pions requires fur- 
ther studies. A possible evidence of the reaction discussed here is a signal from both ZDC's 
and no signal in the central detector. 

In Figj6] we present rapidity distributions of pions y-„+ and rapidity distributions of neu- 
trons t/n- Please note a very limited range of rapidities shown in the figure. The contributions 
for individual diagrams a) - i) (see Fig{T]) are also shown. The diagram d) (from Fig{T]) gives 
the largest contribution. One can observe specific symmetries between different contribu- 
tions on the left and right panels. For instance the long-dash-dotted line on the left panel 
(corresponding to diagram b) ) is symmetric to the dashed line on the right panel (corre- 
sponding to diagram c) ). Clearly, a significant interference effect can be seen. There is no 
region of either pion or neutron rapidity where the diagram (a) dominates. This makes the 
possibility of extracting of tt+tt"'" elastic scattering very difficult. 

For completeness in FigJT] we show the contribution of the diagrams with subleading 
charged reggeon exchanges (see Fig 12]) which could not be seen in the previous plot. We 
show results for the RHIC (left panel) and LHC (right panel) energies. In contrast to the 
other mechanisms the corresponding contribution is rather flat over broad range of rapidities. 
The cross section corresponding to this mechanism is bigger by 2 orders of magnitude for 
the RHIC energy compared to the LHC energy, but rather small compared to the dominant 
contributions shown in FiglH In addition we show contribution of diagrams of FiglS] They 
are comparable to those of diagrams shown in Figj2] at midrapidities but much smaller than 
those from Fig{T] at larger rapidities. We show results of diagrams from Fig{T] with different 
values of the form factor parameter A = 0.8 GeV (bottom dashed line) and A = 1.2 GeV 
(upper dashed line) in order to demonstrate the cross section uncertainties. 

In FigJHlwe present rapidity distributions of pions y.„ for double charged reggeon exchanges 
at W = 500 GeV (left panel) and W = 7 TeV (right panel). The bold solid line represent 
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FIG. 5: Differential cross section da/dr] for neutrons (solid lines) and pions (dotted lines) at the 
center-of-mass energies W = 0.5, 0.9, 2.36, 7 TeV. The smaller bumps include absorption effects 
calculated in a way described in subsection III CI In this calculation we have used A = Aoff = 1 
GeV. The vertical lines at rj = ±8.5 are the lower limits of the CMS ZDC's. The details about 
RHIC ZDC's can be found in Ref.[25]. 

the coherent sum of all amplitudes corresponding to diagrams in Figj2l The contributions 
for individual diagrams are also shown separately. The diagram a) in Figj2] gives the largest 
contribution (long-dashed line). The a2 — F — a2 exchange corresponds to the long-dashed- 
dotted line. One can see that the double reggeon exchange mechanisms shown in Figj2] 
populate midrapidities of the pions and therefore can be measured either at RHIC or at 
LHC In Table HVl we have collected cross section for this component separately for double 
spin conserving (DSC), single spin flip (SSF) and double spin flip (DSF) contributions. All 
this spin contributions are of similar size. The total contribution is about half of nb at RHIC 
(500 GeV) and a few pb at LHC (7 TeV). 

Can the much smaller contribution of diagrams with subleading charged reggeon ex- 
changes be identified experimentally? In Figj9] we show two-dimensional distribution in 
(l/3) 1/4) space. The double-charged reggeon-exchange components from Figj2]are placed along 
the diagonal y^ = 1/4 while the other contributions some distance from the diagonal. There- 
fore imposing 2-dim cuts in the (1/3,1/4) space one could separate the small double charged 
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FIG. 6: Differential cross sections da/dy^+ and da/di/n at W = 7 TeV. Tlie bold solid line rep- 
resent the coherent sum of all amplitudes. The long-dashed (black), long-dash-dotted, dashed, 
long-dashed (red online), dash-dot-dot-doted, dotted, dash-dotted, dash-dot-doted, long-dashed 
(blue online) lines correspond to contributions from a) - i) diagrams. The red, black and blue 
lines correspond to diagrams when neutron, pion and proton are off-mass-shell, respectively. No 
absorption effects were included here. 



reggeons contribution. A very good one- dimensional observable which can be used for the 
separation of the processes under disscusion could be differential cross section da/di/diff, 
where ydiff = Us — Vi and experimentally charged pions should be taken at random (see 
FigHni y IT, first = 1/3 or 1/4 and yn,second = 1/4 or 1/3). For comparison we show contribution of 
diagrams shown in Fig 131 

In FigUll we show distribution of neutrons and pions in the Feynman variable xp = 
2p\\/ ^/s. In this observable the neutrons and pions are well separated. The position of peaks 
is almost independent of energy. While pions are produce at relatively small xp the neutrons 
carry large fractions of the parent protons. The situation is qualitatively the same for all 
energies. 

The distribution in pion-pion invariant mass is shown in FiglT2l Unique for this reaction, 
very large two-pion invariant masses are produced (see e.g. Ref.[16|). The larger energy 
the larger two-pion invariant masses (left panel). The absorption effects almost uniformly 
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FIG. 7: Differential cross sections da/dy^+ sX W = 500 GeV (left) and M^ = 7 TeV (right). The 
lines represent the coherent sum of all amplitudes from diagrams in Fig(T]-[3l as well as separate 
classes of contributions marked by the number of the figure where they are shown. No absorption 
effects were included here. 

TABLE IV: Full-phase-space integrated cross section (in nb) for exclusive 7r"'"7r"'" production for the 
amplitude with the double charged reggeon exchanges (diagrams in Figj2]) at the center-of-mass 
energies W = 0.5, 7 TeV. No absorption effects were included here. The meaning of the acronyms: 
DSC - double spin conserving, SSF - single spin flip, DSF - double spin flip. 



exchange 


W = 0.5 TeV 


W = 7 TeV 


p+ _ ^0 _ p+ 


0.43 


3.3x10-3 


P^-al-p+ 


0.14 


1.0x10-3 


at — p^ — at 


0.11 


5.4x10-4 


p~^ — LO — p'^ 


1.5x10"^ 


1.1x10-*^ 


sum of all amplitudes 


0.7 


5.1x10-3 


DSC 


0.17 


1.5x10-3 


SSF 


0.18 


1.3x10-3 


DSF 


0.18 


1.0x10-3 


at -IP -at 


4.4x10-3 


2.5x10-3 



reduce the cross section. We show also distributions with different values of the form factor 
parameter in order to demonstrate the cross section uncertainties (righ panel). 

The distributions in the transverse momentum of neutrons and pions are shown in FigJTSl 
The figure shows that the typical transverse momenta are rather small but large enough to 
be measured. The distributions for neutrons are rather similar to those for pions. 

The energy distributions of neutrons are presented in FigJT4l Generally the larger colli- 
sion energy the larger energy of outgoing neutrons. When combined with the previous plot 
it becomes clear that the neutrons are produced at very small polar angles (large pseudora- 
pidities) and can be measured by the ZDC's (see also Figj5]). 
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FIG. 8: Differential cross sections da /dy-,^ for double charged reggeon exchanges a,t W = 500 GeV 
(left) and W = 7 TeV (right). The contributions for individual diagrams in Figl2] are shown 
separately. No absorption effects were included here. 
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FIG. 9: Differential cross sections in (2/3,2/4) space at W = 500 GeV (left) and W = 7 TeV (right). 
The coherent sum of all amplitudes from diagrams in FiglH Figl3]and the contribution of diagrams 
in Figl2] with double-exchange reggeons placed along the diagonal are presented. No absorption 
effects were included here. 



In FiglTS] we show two-dimensional correlations between energies of both neutrons mea- 
sured in both ZDC's. The figure shows that the energies of both neutrons are almost not 
correlated i.e. the shape (not the normalization) of da/dE^ {da/dEn2) is almost indepen- 
dent of i?„2 (Eni)- There should be no problem in measuring energy spectra of neutrons on 
both sides as well as two-dimensional correlations in {En^,En2)- 

Finally in FigJT6]we present the distributions in azimuthal angle (f) between the transverse 
momenta of the outgoing neutrons (pions) . Clearly a preference of back-to-back emissions 
can be seen. The measurement of azimuthal correlations of neutrons will be not easy with 
first version of ZDC's as only horizontal position can be measured. Still corellations of 
horizontal hit positions on both sides could be interesting. A new correlation observable, 
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FIG. 10: Differential cross sections da /dy diff at W = 500 GeV (left) and VF = 7 TeV (right). The 
lines represent the coherent sum of all amplitudes from diagrams in FigtTJ Figl3]and the contribution 
of diagrams in Figl2]with double-exchange reggeons placed at ydiff ~ 0. No absorption effects were 
included here. 



tt. 5 
X 

-a 

^ 4 



W = 7TeV, 





-0.5 



0.5 



FIG. 11: Differential cross section da/dxp for the pp 
absorption effects were included here. 



nmr^TT^ reaction at W = 7 TeV. No 



taking into account possibilities of the apparatus, should be proposed. In contrast the two 
7r"'"'s are almost not correlated in azimuthal angle. However, such a distribution may be not 
easy to measure. 

We have shown that at present the reaction under consideration can be strictly measured 
only in a rather limited part of the phase space (midrapidities of pions) where the cross 
section is rather small and where the double charged reggeon mechanism dominates. In 
Table |V] we have collected the cross sections in nb for different experiments at LHC and 
RHIC At LHC where the separation of the double-reggeon exchange mechanism is possible 
the cross section is rather small of the order of a fraction of pb. At RHIC the cross section 
with experimental cuts should be easily measurable as it is of the order of a fraction of nb. 
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FIG. 13: Differential cross section da/dpt for the pp — )• nnvr^vr^ reaction at W = 7 TeV. The solid 
and dotted lines correspond to the distribution in the transverse momentum of neutrons and pions, 
respectively. The lower curves correspond to calculations with absorption effects. 



IV. CONCLUSIONS 

We have estimated cross sections and calculated several differential observables for the 
exclusive pp — > nmT~^7c^ reaction. Because our parameters are extracted from the analysis 
of known two-body reactions we expect that our prodictions of the cross section are fairly 
precise inspite of the complications of the reaction mechanism. The full amplitude was 
parametrized in terms of leading pomeron and subleading reggeon trajectories. We have 
consider 3 classes of diagrams. The first class gives the largest contribution but concen- 
trated at forward or backward pion directions. There are also diagrams with double charged 
exchanges with subleading reggeons p"*" and a^. Although the cross section for these contri- 
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butions is rather small, it is concentrated at midrapidities of pions where the cross section 
can be easily measured. The double-exchange reggeons processes can be separated out in the 
two-dimensional space of rapidities of both pions or in the distribution of the pion rapidity 
difference. 

Large cross sections have been obtained, even bigger than for the pp — )■ pp-K^-K^ reaction 

Iq . Several mechanisms contribute to the cross section, which leads to an enhancement 

of the cross section due to interference effects. These interference effects cause that the 

extraction of the elastic Tr+Tr"*" cross section as proposed recently [31] seems in practice rather 

impossible. 

The specifity of the reaction is that both neutrons and pions are emitted in very for- 
ward/backward directions, producing a huge rapidity gap at midrapidities. While the neu- 
trons could be measured by the ZDC's, the identification of pions may be difficult. We 
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FIG. 16: Azimuthal angle correlations between neutrons and beetwen pions for the pp — >• nmr~^TT^ 
reaction at W = 7 TeV. The lower curves correspond to calculations with absorption effects in- 
cluded. 

TABLE V: Cross section (no absorption effects) with different experimental cuts on pt^-K, ^tt and 





W (TeV) 


Pt,7T > 


IvttI < 


IVnlzDC ^ 


a (nb) 


ALICE 


7 


0.15 


0.9 


8.7 


6.3x10^5 


ALICE 


7 


0.15 


1.2 


8.7 


1.2x10-^ 


ATLAS 


7 


0.5 


2.5 


8.3 


4.9x10"^ 


CMS 


7 


0.75 


2.4 


8.5 


4.5x10"^ 


RHIC 


0.5 


0.2 


1 


- 


2.0x10-2 



think that the measurement of both neutrons and observation of large rapidity gap is a very 
good signature of the considered reaction. We expect the cross section for the nnTr+Tr+Tr", 
nnn^Ti^Ti^Ti^, etc., which could destroy rapidity gaps, to be smaller but a relevant estimates 
need to be done. In addition for events with larger number of pions the rapidity gap would 
be destroyed. Therefore the formally kinematically incomplete measurement of two neutrons 
only could be relatively precise. We have found that the neutrons measured in ZDC's seem 
to be almost uncorrelated in energies. 

We have made predictions for azimuthal angle correlations of outgoing neutrons. Such 
distribution should be possible to measure in a future. At present at CMS only horizontal 
position can be measured. We have predicted back-to-back correlations with a sizeable 
diffusion. 

We have included elastic rescattering effects in a way used recently for the three body 
processes. These effects lead to a substantial damping of the cross section. The bigger 
energy the larger the effect of damping. Other processes (e.g. inelastic intermediate states 
or final state vr+ra interactions) could lead to additional damping. At present there is no full 
understanding of the absorption effects. A future experiment could provide new data to be 
analyzed and could shed new light on absorption effects which are essential for understanding 
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exclusive processes, even such important ones as exclusive production of the Higgs boson. 

In the light of our analysis it becomes clear that extraction of the elastic 7r"'"7r"'" cross 
section seems impossible, due to interference of several processes discussed in our paper. We 
did not find any corner of the phase space where the relevant diagram dominates. 

There is an attempt to install forward shower counters in the LHC tunnel. Most probably 
they will not be able to measure energy of the pions but they can signal some activity there. 
We expect that "some activity" will mean, with a high probability, just one vr"*" on one side 
and the other vr"*" on the other side. 



V. APPENDIX 
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FIG. 17: Diagrams for various exchanges in irp collisions. 

The p-meson/reggeon and a2-meson/reggeon exchanges are known to have not only the 
nucleon spin-conserving part but also the dominant nucleon spin-flip component while the 
cu-meson/reggeon exchange to nucleons is mainly spin-conserving. We write the amplitude 
for the reggeon exchanges (see Fig lTTj) in the following compact phenomenological form: ^ 



jreggeon—exch 



MlT:T-Z-is,t) 



\ Qi(t)-1 / p 

S \ I iJMN 



X r]isCl[ — ] exp ( -^^{t - t^i,^) ) 5\x,j\i , (5.1) 



and the pion exchange amplitude as 

-^C^lt, A^l-^'^) = 9-kNN F^NNit) u{Pn', Xn'YiiMpn, Atv) 



i /s\""^*^ 

X (A;^ + qne;{k', A,) g,^^ F,^^{t) - . (5.2) 



Above the \/ —{t — tmin)/Mo factor is due to the meson spin-flip (in the tt — >■ cj, tt — >■ p 
and TT — )■ 02 transitions), Mq is a reference scale factor taken here Mq = 1 GeV (which is used 
here to have the same units for the coupling constants). The double spin-flip components 
do not interfere with the spin-conserving ones and can be calculated separately. Here we 
have introduced one more phenomenological (dimensionless) parameter r^ which describes 
coupling for the spin-flip components. It is known to be of r^ = 7.5, r^ ~ 6.14, r^ ~ 0.17 
|22| and r^ ~ 8, r^ ~ 4.7, r^ ~ 0.9 |23|. In the present calculations we take r^ = 7.5, r^ 



^ For the case of the tt p — > uj^n, tt p — > p'^n and tt p — )• CjTi reactions the amplitude should be multiplied 
by \/2 which is relatad to isospin Clebsch-Gordon coefficient. 
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= 6 and r^ = 0. The coupling constant gp.,,-,, is taken as g'^^^/iTr = 2.6. The form factors 
are parametrized as F{t) = exp {{t — m'^) / A"^) . We improve the parameterization of the 
amphtude fl5.2p by multiplying by the factor (s/sq)""^*-', where aT,{t) = a'^{t — m^) is the 
pion Regge trajectory with the slope of trajectory a'^ = 1 GeV~^. 

We adjust the C[ (where i = ]P,p,ijj,a2) coupling constants to the world experimental 
data often obtained from partial wave analysis in the three-pion system. The effective 
normalization constants for the auxiliary reactions are related to those in the A^A^ scattering 
and the g],^a2,p,oj coupling constants we need in our problem as: 



CI = sjcr ■ gU, . (5.3) 

Since Cf^ are known from phenomenology (Table [T|), 5'^_j.j can be obtained from our fits: 
^f_,,^ = 1.4 GeY-\ gP^^,^ = C_, = 22 GeV-^ and g^^^ = g^^^ = 4 GeV^ 

In Fig US] we show the total cross section for the 7i~p — )■ Oj p, tt~p — )■ u^n, n~p — )■ p^n and 
TT^p — )■ p^p reactions as a function of the incident-beam momenta Piab- Our fit is shown by 
the solid line. In the panel a) {ir^p — )■ a2P reaction) we show individual contributions of p 
and pomeron exchanges. The pomeron exchange dominates at high energies whereas the p 
exchange at small energies. This separation of mechanisms allows to extract two independent 
coupling constants. We show also spin-conserving and spin-flip amplitudes separately. In 
panel b) we show our fit for the vr'p — )■ u^n. Here only p exchange is possible. In panel 
c) (tt^p — i- p^n reaction) we show contributions for charged pion exchange (parameters 
fixed from phenomenology) and 02 exchange (parameters found from the analysis of the 
n~p — )■ a^p (see panel a))). Finally in panel d) (vr^p — )■ p'^p reactions) we show contributions 
for neutral pion exchange, 02 exchange and u exchange (relevant coupling constant found 
from the analysis of the n^p — )■ u^n reaction (see panel b))). 

Having fixed the parameters we can proceed to our four-body pp — t- nmi^Ti^ reaction. 
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